Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 5 February 2008 (MN KTeX style file v2.2) 



SN 2005CS in M51 

I. The first month of evolution of a subluminous SN II plateau 

. . . A. Pastorello^ *, D. Sauer^, S. Taubenberger\ P. A. Mazzali^'^'^'"^, K. Nomoto^ '^, 
O K. S. Kawabata^, S. Benetti^, N. Elias-Rosa^, A. Harutyunyan^, H. Navasardyan^, 
<N L. Zampieri^ T. lijima^ M. T. Botticella^ ^ G. Di Rico^ ^ M. Del Principe^ 
^ M. Dolci^, S. Gagliardi^, M. Ragni^'^, G. Valentini^ 

' Max-Planck-Institut fiir Astrophysik, Karl-Schwarzschild-Str. 1, 85741 Garching bei Munchen, Germany 
^— ( , - INAF Osservatorio Astronomico di Trieste, Via Tiepolo 11, 34131 Trieste, Italy 

■ ^ Department of Astronomy, School of Science, University of Tokyo, Bunkyo-ku, Tokyo 113-0033, Japan 
' Research Center for the Early Universe, School of Science, University of Tokyo, Bunkyo-ku, Tokyo 113-0033, Japan 
' Hiroshima Astrophysical Science Center, Hiroshima University, Hiroshima 739-8526, Japan 
I , * INAF Osservatorio Astronomico di Padova, Vicolo dell' Osservatorio 5, 35122 Padova, Italy 
^ : ' INAF Osservatorio Astronomico di Padova, Sezione di Asiago, Via dell'Ossen'atorio 8, 36012 Asiago (Vicenza), Italy 
Universita degli Studi di Teramo, Viale Crucioli 122, 64100 Teramo, Italy 
^— ^ ' INAF Osservatorio Astronomico di Collurania, via M. Maggini, 64100 Teramo, Italy 

■ 

in ; 

C ^ ^ Accepted ; Received in original form .... 

o 



O 



ABSTRACT 

Early time optical observations of supernova (SN) 2005cs in the Whirlpool Galaxy (M51), 
are reported. Photometric data suggest that SN 2005cs is a moderately under-luminous Type 
II plateau supernova (SN IIP). The SN was unusually blue at early epochs (U-B ^ -0.9 about 
three days after explosion) which indicates very high continuum temperatures. The spectra 
show relatively naiTow P-Cygni features, suggesting ejecta velocities lower than observed in 
more typical SNe IIP. The earliest spectra show weak absorption features in the blue wing of 
the He I 5876A absorption component and, less clearly, of H/3 and Ha. Based on spectral 
modelling, two different interpretations can be proposed: these features may either be due to 
high-velocity H and He I components, or (more likely) be produced by different ions (N II, Si 
II). Analogies with the low-luminosity, ^''Ni-poor, low-velocity SNe IIP are also discussed. 
While a more extended spectral coverage is necessary in order to determine accurately the 
properties of the progenitor star, published estimates of the progenitor mass seem not to be 
consistent with stellar evolution models. 

Key words: supernovae: general - supernovae: individual (SN 2005cs) - supernovae: individ- 
ual (SN 1997D) - supernovae: individual (SN 1999br) - supernovae: individual (SN 2003Z) - 
galaxies: individual (M51) 



1 INTRODUCTION 

Type II supernovae (SNe II) are believed to be produced by the ex- 
plosion following the core-collapse of massive stars that retained 
most of their H envelope at the time of explosion. Some SNe II 
spend a period at almost constant luminosity: this phase, lasting 
sometimes a few months, is called "plateau" (hence the label SN 
IIP). When the SN enters this phase, the temperature is low enough 
that the massive H envelope, initially ionized because of the de- 
position of energy by the shock wave, starts to recombine. After 
recombination, the light curve of SNe IIP declines steeply, until 
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it settles onto the "radioactive tail", as do other SN Types. In this 
phase the luminosity is due mainly to the radioactive decay of ^'Ni 
to ^''Co to ^''Fe. 

Despite t he considerable number of SN e IIP studied in re- 
cent years te.g. lPatat et alll9'94t ^Hamuv 2003), the physical prop- 
erties of the progenitor stars are still a matter of debate. Thanks 
to the direct identification of several SN precursors in deep pre- 
explosion images, important constraints have been established on 
the nature of the progenitors of SNe IIP. The first SN with a 
known progenitor was SN 1987A in the Large Magellanic Cloud. 
Archival images showed th at the sup ergiant progenitor was unusu- 
ally blue (e.g. Sonnebom et al .119871) . The present e nsemble of SNe 
IIP with detected progenitor includes SN 2003gd tVan Dvk et alJ 
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Q^zAMS = 9*1 Mb) of spe ctral type in the range K0-M4 \Ia et alj 
l200d iMaund et alj|200 j) . These two papers report different val- 
ues for the I band detection magnitude: 24.15 (i.e. absolute magni- 
tude M; » -5.5) an d 23.3 (M; x - 6.4, adopting the same redden- 
ing and distance as lLi et all i200^). r espect ively. The position of 
the candidate as measured bv lLi etalJ <200dl is a = 13''29'"52!76, 
6 = +47°10'36'.'ll (J2000). Alternatively, Richmond (2005) claims 
that SN 2005cs exploded near a cluster of young stars and finds that 
the progenitor candidate could be a blue star at a = 13''29"'52'!803, 
S = +47°10'36'.'52 (J2000), with Mv « -6. However, since 
SN 2005cs is evolving like a normal plateau event (see Sect. |2j, 
this blue supergiant progenitor candidate is not convincing. 

Here we present early time optical observations of SN 2005cs 
(until ~ 1 month after the explosion). In Sect. 2 we describe the 
photometric evolution of SN 2005cs and in Sect. 3 we analyse the 
spectroscopic data. A discussion follows in Sect. 4. 




Figure 1. SN 2005cs in M51: V band image obtained on July 1, 2005 with 
the Copernico 1.82m t elescope o f Mt. Ekai\ Asiago (Italy). The sequence 
stars from^Richmond et alj ll99d) are labelled with letters. Some of our lo- 
cal sequence standards (Pastorello et al. in prep.) are indicated by numbers. 



20031 : ISmartt et all2004lHendrv et alliooj) . SN 2004et <Li et alj 



2005ft. SN 199 9ev (Maund & Smartt 2005) and, now, SN 2005cs 
Li et al-20 



iMaund et all20Q5:,Li et al-2006.) . Lower magnitude limits or am- 
biguous detections ha ve been obtained for ot her SNe IIP precur- 
sors (e.g. SN 1999br. iMaund & Smartt"2005'). All these observa- 
tions seem to support the idea that most SNe IIP originate from the 
explosion of moderately massive stars (M < 15o). 

SN 2005cs was discovered in the famous Whirlpool Galaxy 
(NG C 5194 orM5 1) by Kloehr et al. (2005) on 2005 June 28.905 
UT. Modiaz et al. (2005) classified it as a young Type 11 SN. 
The earliest detection was that of M. Fiedler on Ju ne 27.91 UT 
(SNW eb'). Nothing was visible on June 20.6 UT jKloehr et all 
l2005h at the SN position. Moreover, no clear detection of SN 
2005cs was found on several different images obtained on June 26 
by other amateur observers (SNWeb). In particular, the SN site was 
monitored by the team of the Osservatorio Astronomico "Gemini- 
ano Montanari" (Cavezzo, Modena, Italy) on June 26.89 using a 
Newton 0.4m telescope and nothing was detected below the follow- 
ing limits: B ^ 17.3, V ^ ill, R ^ 17.6. These detection limits con- 
strain the explosion time to a very small uncertainty (about 1 day). 
Therefore in this paper we adopt June 27.5 UT (JD = 2453549+1) 
as the explosion epoch . 

The coordinates of SN 2005cs are a = 13''29"'5Z85 and 
d = -i-47°10'36'.'3 (J2000). The object lies in the southern arm 
of M51, 15" West and 67". 3 South of the galaxy nucleus (Fig. 
[T}. M51 is classified by NED^ as an SA(s)bc peculiar galaxy. 
The galaxy also hosted the well studied core-collapse SN 19941 
fWheeler et al.'l994V'Fili Dpenko et alJl995tlRichmond et aljl996t 
jciocchiatti et al. 1996). 

A candidate progenitor for SN 2005cs was identified in com- 
bined HST ACS F814W images (~ 1 band) as a red supergiant 



2 PHOTOMETRY 

Our photometric data were obtained using seven difi"erent tele- 
scopes, and cover 19 epochs (including the prediscovery limit), un- 
til approximately 35 days after the explosion. 
All data were pre-reduced with standard IRAF^^ procedures, and in- 
strumental SN magnitudes were determined using the point-spread 
function (PSF) fitting technique performed with the "SNOOPY"'* 
package. Since SN 2005cs is a bright object, this technique pro- 
vides acceptable results, although the background region of SN 
2005cs is extremely complicated, such that the subtraction of the 
template could be more appropriate when the SN fades. 

In order to transform instrumental magnitudes to the standard 
Johnson-Cousins system, first-order colour corrections were ap- 
plied with colour ter ms deriv ed from observations of photometric 
standard fields ( Landoltll992h . The photometric zeropoints were fi- 
nally determined by comparing the magnitudes of a local sequence 
of stars in the vicinity of M51 (cf. Fig. to the values reported 
for some of these stars by Richmond et al. (1996) in their study on 
SN 19941. A complete definitive photometry (calibrated on a larger 
sequence of stars, i.e. those labelled by numbers in Fig.0 will be 
presented in a forthcoming paper (Pastorello et al. in prep.). 

The SN magnitudes are reported in Tab. Q and the U, B, V, 
R, 1 light curves are shown in Fig.Otogeth er with th e light curves 
of SNe 1999br jHamuvll200ltlPa storello e t alj|2004 and 1999em 
faamuv et all 1200 iL Leonard et alj i2002l: Elmhamdi et"ai]|2003h . 
shifted arbitrarily in magnitudes (but not in phase) in order to 
match the light curves of SN 2005cs. The overlap of the B, V, R 
and 1 band light curves of the three SNe is very good, while some 
differences are visible in the U band evolution: the U band light 
curve of SN 2005cs declines more rapidly than that of SN 1999em. 



http://www.astmsurf.comJsnweb2/\ 



http://nedwww. ipac. caltech. edu/\ 



According to lModiaz et alj J2005h . the simultaneous presence 
in the SN spectrum of narrow Galactic and host galaxy interstel- 
lar Na ID lines with analogous equivalent width (0.2 A) indicates 
a similar contribution of the two components to the total extinc- 
ti on. The host galaxy e xtinction can be estimated us ing the relation 
of lTuratto et alj |20o3), while the value provided bv lSchlegel et alJ 



IRAF is distributed by the National Optical Astronomy Observatories, 
which are operated by the Association of Universities for Research in As- 
tronomy, Inc, under contract to the National Science Foundation. 

SNOOPY is a package implemented in IRAF by E. Cappellaro, based on 
DAOPHOT. 
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Table 1. UBVRI magnitudes of SN 2005cs and assigned errors. Both measurement errors and uncertainties in tlie pliotometric calibration (tlie most important 
source of errors) are taken into account. Measurement uncertainties give a minor contribution to the total error because, despite the relatively complex 
background, the SN was much brighter than any other nearby source. 



Date 


JD 

(+2400000) 


U 


B 


V 


R 


/ 


Inst. 




S^54.S 




> 17 3 


17 7 


^ 17 6 




NCO 


ns/n(S/'?n 


53557 3fi 


1 3 48+0 05 


14 3ft+n 05 




1 4 4fi+0 04 


1 4 44+0 04 


Caha 


VJJ/ W / /U i 


J J J J J . -J J 


1 1. S^+0 0^1 


]A -^A+O OS 


14 4(S+0 O'S 


1 4 40+0 0? 


1 4 ^4+0 Ofi 




05/07/09 






14. 4.S+0 04 


14 S 1 +0 0^ 


14 47+0 0^ 


14 9q+n 10 


Ekar 


05/07/05 


53557.42 




14.52+0.04 


14.54+0.04 


14.42+0.05 


14.29+0.08 


Ekar 


05/07/06 


53557.84 


14.02+0.06 


14.60+0.05 


14.55+0.03 


14.36+0.03 


14.40+0.03 


Sub 


05/07/07 


53559.40 




14.65+0.04 


14.56+0.09 


14.41+0.09 


14.35+0.11 


TNT 


05/07/1 1 


53563.38 


14.81+0.05 


14.87+0.04 


14.56+0.04 


14.37+0.02 


14.23+0.05 


Ekai' 


05/07/11 


53563.42 


14.89+0.04 


14.93+0.04 


14.53+0.04 


14.37+0.03 


14.26+0.03 


TNG 


05/07/13 


53565.38 


15.27+0.06 


15.04+0.05 


14.64+0.02 


14.37+0.03 


14.29+0.05 


LT 


05/07/14 


53566.36 


15.29+0.05 


15.11+0.05 


14.68+0.02 


14.44+0.04 


14.23+0.02 


Ekar 


05/07/14 


53566.40 




15.09+0.04 


14.66+0.04 


14.40+0.05 


14.32+0.06 


TNT 


05/07/17 


53569.42 


15.92+0.07 


15.31+0.05 


14.67+0.03 


14.40+0.03 


14.27+0.04 


LT 


05/07/19 


53571.40 




15.36+0.12 


14.72+0.05 


14.45+0.09 


14.26+0.07 


TNT 


05/07/20 


53572.40 




15.39+0.07 


14.71+0.03 


14.45+0.03 


14.27+0.03 


TNT 


05/07/25 


53577.40 




15.46+0.09 


14.72+0.04 


14.36+0.06 


14.24+0.06 


TNT 


05/07/27 


53579.40 




15.60+0.07 


14.73+0.04 


14.39+0.04 


14.16+0.04 


TNT 


05/07/31 


53583.39 


16.83+0.09 


15.78+0.06 


14.71+0.03 


14.38+0.03 


14.12+0.05 


TNG 


05/07/31 


53583.47 


16.81+0.11 


15.77+0.06 


14.69+0.03 


14.36+0.03 


14.15+0.04 


LT 



NCO = Osservatorio Cavezzo 40cm Newton Telescope 

Caha = Calar Alto 2.2m Telescope + CAFOS 
Ekar = Asiago 1.82m Copemico Telescope + AFOSC 

Sub = Subaru Telescope 8.2m + FOCAS 
TNG = Telescopic Nazionale Galileo 3.5m + Dolores 
LT = Liverpool Telescope 2.0m + RATCAM 
TNT = Teramo Normale Telescope 72cm 

il998h is adopted for the Galactic extinction. Taking into account 
both components, a total reddening of E(B-V) = 0.06 is estimated. 
However, other methods yield slightly higher values for the red- 
dening. The study of the nearby H II region CCM 56 and the 
three-colour photometry of some red supergiants close to the SN 
positi on yield a total redd ening of E(B-V) = 0.16 Jsresolin et al] 
|2004^ and E(B-V) = 0.12 iMaund et all2005l) . respectively Lack- 
ing strong constraints in favour of one of these three methods, an 
average value of E(B-V) = 0.11 ±0.04 is adopted. However, this 
relatively low amount of interstellar reddening is consistent with 
the blue colour of the early spectra of SN 2005cs. 

iFeldmeier et all il997h used planetary nebulae to determine a 
distance to M5 1 of 8.4 Mpc (i.e. distance modulus yu = 29.62±0. 15). 
Averaging all R band data reported in Tab.Q an R band absolute 
magnitude during the plateau phase Mjt x -15.48±0.16 was ob- 
tained. This implies that SN 2005cs was relatively underluminous 
compared to more t ypical SNe IIP and w as therefore similar to SNe 
1997D and 1999br iTuratto et al.'l998':'Benetti et al.''200l':l Hamuvl 
I2OOI: Zampieri et al..2003.: Pastorel lo et al. 2004), although less 
extreme. 

In Fig. we compare the early time B-V (top) and V-R 
(middle) colour curves, and the pseudo-bolometric uvoir light 
curve (bottom) of SN 2005cs with those of t he SNe IIP 1999br 
('Ha muvl200i[ |Pastorello et al.'2004'). 1999em faamuv et all200lt 
Leon ard etall20 02: Elmhamdi et al. 2003) and the peculiar 1987A 
(Men zies et all 11987). The light curve of another moderately 
under-luminous object, SN 2002gd I Pastorella.2003k) . is also in- 
cluded in Fig.0bottom) for comparison. 

We adopted a distance modulus yu = 30.97 and a total B band 
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Figure 2. Early U, B, V, R, I light curves of SN 2005cs, including our limits 
of June 26 (when the SN was not detected) and the unfiltered limit from lAU 
Circ. 8853 (asterisk). Also the light curves of SN 1999br (BVRI) and SN 
1999em (UBVRI) have been included for comparison, and shifted in mag- 
nitude by an arbitrary amount in order to match the corresponding curves of 
SN 2005cs. The very eariy B band detections of M. Fiedler (SNWeb) have 
been also reported (filled pentagons). 
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Figure 3. Top: early B-V colour curve of SN 2005cs and those of SNe 1999em, 1999br and 1987A. Middle: same as above, but for V-R colour Bottom: early 
uvoir pseudo-bolometric light curves for the same sample of SNe IIP and, in addition, SN 2002gd. For references, see the text. 



extinction Ah 



= 0.1 (Pastorello et al."2004') for SN 1999br, ^i = 
30.34 jLeonard et a"h,.2003) and Ag = .41 (Baron et al.,2000:.) for 



SN 1999em, n = 18.49 and Ag = 0.79 iAmett et al Il989t) for SN 
1987 A, and ^ = 32.87 and Ag = 0.29 |Pastorelldl2003h for SN 
2002gd. Distance moduli for SNe 1999br and 2002gd were esti- 
mated from the recession velocity corrected for the Local Group 
infall onto the Virgo Cluster (vv,, ), adopting a value of Ho = 72 km 
s"' Mpc"'. 

All typical SNe IIP show a similar colour evolution. Their 
colour curves become monotonically redder with time. On the other 
hand, SN 1987 A was significantly redder at early epochs (until ~ 
50 days), but then the colour difference between it and other SNe 
IIP decreases with time. SN 2005cs follows the behaviour of SNe 
IIP: the B-V colour increases from -0.2 to 1 during the first month, 
while the V-R colour increases much more slowly (0 to 0.3) over 
the same period (Fig.|3|top and middle). 

Interestingly, the uvoir light curve of SN 2005cs has an 
evolution reminiscent of that of normal SNe IIP (Fig. |3| bot- 



tom). The uvoir luminosity of SN 2005cs is intermediate between 
those of the interme diate ^*Ni mass SN 1999em (M^, ~ 0.05Mo, 
Leonard et all2003h and the low-luminosi t y, ^''Ni-poor SN 1999br 
(Mm « 2xlO~^M^. Fzampieri et ai]l2003t1Pastorello et alj|2004b . 
and is similar to that of the '*Ni-poor SN 2002gd (see Sect. Hand 
|Pastorellcl200l . 



3 SPECTROSCOPY 

Spectroscopy is available for 10 different epochs, ranging from 
about 3 to 35 days after the explosion. A summary of all spectro- 
scopic observations is given in Tab.|2| 

All raw frames were first bias and flat-field corrected, and 
then the SN spectra were optimally extracted. Wavelength calibra- 
tion was obtained with the help of comparison lamp exposures, 
while the spectra were flux calibrated using standard star spec- 
tra obtained on the same night as the SN observations. When no 
spectro-photometric standard star was observed, a sensitivity func- 
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Table 2. Journal of spectroscopic observations of SN 2005cs. The phase is from the explosion epoch. 



Date 


JD 
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Resolution" 




+2400000 
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(A) 


05/06/30 
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3 
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9 


05/07/01 
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4 
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05/07/02 
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05/07/05 
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8 
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9 
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14 
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3500-7800 


24 


05/07/11 


53563.44 


14 


TNG+DOLORES+gm.LRB,LRR 


3150-9500 


18,17 


05/07/14 


53566.35 


17 


Ekar+AFOSC+gm.4,2 


3500-9200 


24,38 


05/07/15 


53567.42 


18 


Pennar+B&C+150tr/mm 


4150-8750 


25 


05/07/19 


53571.45 


22 


Pennar+B&C+ 1 50tr/mm 


4050-7900 


25 


05/07/31 


53583.44 


34 


TNG+DOLORES+gm.LRB,LRR 


3150-9600 


18,17 



" as measured from the full-width at half maximum (FWHM) of the night-sky lines 
CA3.5m = 3.5m Telescope, Calar Alto Obs., Centro Astr. Hispano Aleman, Almeria (Spain) 
Ekar = 1.82m Copemico Telescope, INAF - Osservatorio di Asiago, Mt. Ekar, Asiago (Italy) 

Pennar = 1.22m Galilei Telescope, Universita di Padova, Loc. Pennar, Asiago (Italy) 
Subaru = 8.2m Subaru Telescope, National Astr Obs. of Japan, Mauna Kea, Hawaii (USA) 
TNG = 3.5m Telescopio Nazionale Galileo, La Palma, Canary Isl. (Spain) 



tion derived on a different night (close in time) was used. Telluric 
features were removed from the SN spectra, again using spectro- 
photometric standard spectra. However, the strong telluric band at 
7570-7750A is superimposed on the O I AlllA SN feature, and im- 
perfect removal might significantly affect the profile of this line. Fi- 
nally all spectra were checked against the quasi-contemporaneous 
photometry and, where discrepancies occurred (especially during 
low-transparency or bad-seeing nights), the photometric data were 
used to derive a scaling factor to apply to the SN spectrum. The rel- 
ative, final flux calibration was reasonable, and the agreement with 
photometry within 10%. 



3.1 Spectral evolution 

We monitored SN 2005cs spectroscopically for about one month. 
There are only two significant observational gaps, between July 6 
and 11 (phase 9-14 days) and between July 19 and 31 (phase 22- 
34). During the last part of the coverage the spectrum of SN 2005cs 
evolved significantly. The spectral sequence is shown in Fig. El 

The first spectra (phase 3-5 days) are characterized by a very 
blue continuum. The most prominent features are the P-Cygni pro- 
files of the H Balmer lines and of He I 5876A. The position of the 
minimum of these lines indicates expansion velocities of the ejecta 
between 5000 and 8000 km s"', significantly lower than th ose typ- 
icall y observed in SNe IIP at a comparable phase (cf. e.g. lHamuvl 
1200k .2003.: ,PastorellQ..,200S) . The most intriguing property is the 
presence of absorption features on the blue side of Hyf?, of He I 
5876A (very prominent) and, though this is clearly visible only in 
more evolved spectra, of Ha. These absorptions can be interpreted 
as either high-velocity (HV) H I and He I features, or as absorption 
lines of other ions (N II and Si II). The possible presence of HV 
lines could be explained either by an unusual density structure of 
the progenitor star, with a relatively dense and He-rich outer shell, 
or by strong mixing and asphericity of the ejecta. As discussed in 
Sect. l3.2l there is evidence that these lines are due to N II and Si II. 

The subsequent two spectra (phase 8-9 days) also show a blue 
continuum, but the features near 4580A and 5580A have com- 
pletely disappeared. The Fe II multiplet 42 lines (^4924, /15018, 



/i5169) begin to be visible to the red side of Hj8. The He I 
5876A line becomes dimmer and a strong O I AlllA line is now 
visible. Moreover, a very prominent absorption feature is now well 
developed at ~6300A, close to the blue wing of the Ha absorption. 
We tentatively identify this line as the Si II 6347, 6371 A doublet 
(the feature that is prominent in the photospheric spectra of SNe la, 
hereafter Si II 6355A), rather than as a HV Ha component (see also 
Sect.|T2|and Sect.lT^. 

The spectra at phases 14-22 days show redder continua and 
lines with deeper P-Cygni profiles. In the region below ~5300A, 
together with the Balmer H lines, we identify several metal lines. 
In addition to strong lines of Fe II, Ti II and Sc II (e.g. the ab- 
sorptions at 5250A and 5470A), Sr II ^4078, AAlld (doublet 1) 
and /i4162, /i4305 (doublet 3) are possibly detected. Consistently 
with the lower effective temperature, the He I line is no longer vis- 
ible and the feature near 5800A is probably due to the increasing 
strength of Na ID. Finally, Ca II H&K and the Ca II IR triplet are 
now among the most prominent features. The components of the 
Ca II IR triplet are not blended, confirming the low ejecta velocity. 

By the times of the last TNG observation, obtained about 34 
days after the explosion, the spectrum had noticeably changed. 
Overall, it looks like a spectrum of a typical SN II during the re- 
combination phase. The continuum is much redder, and the flux 
deficit below ~3800A is due to strong line blanketing from the Fe 
II features I Mazzali et al. 1992), but other metal ions may also con- 
tribute significantly to the continuum shape. The most prominent 
lines are now Ho', Ca II H&K and the Ca II IR triplet, all with very 
well developed P-Cygni profiles. A detailed line identification of 
this spectrum is presented in Sect. l3.3l 



3.2 High-velocity H I, He I features or N H, Si H lines? 

Fig. |5| shows the early evolution of the spectral region between 
4100A and 7300A. The Hj8, He I 5876A and Ha absorption min- 
ima are marked and the corresponding line velocities are reported. 
The features at 4580A, 5580A and near 6300A are also marked, 
with labelled the two alternative identifications: as putative HV 
H and He I (with the corresponding line velocities), and as N 
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Figure 4. Spectroscopic evolution of SN 2005cs from ~3 days to about one montli after explosion. Telluric features were not removed from the noisy B&C 
spectrum of July 19. For details, see also Tab.l2l All spectra are at the host galaxy rest frame. 



II and Si II lines (bottom of Fig. |5j- iBaron et'al] i200(t) discuss 
the identification of the 4580A and 5580A features in an early 
time spectrum of SN 1999em. Using the simple parametrised code 
SYNOW I Fisher 2000), in which the relative line strengths for each 
ion are fixe d assuming local thermodynamic equilibrium (LTE), 
iBaron et'al! jipOO) find that these features could be consistent 
with N II ^4623 and A5619. However, the non-LTE model at- 
mosphere code PHOENIX provides a synthetic spectrum which 
leads them to reject this identification, because it would require 
an overabundance of nitrogen to reproduce these features. There- 
fore they support the identification of the 4580A, 5580A absorp- 
tions as secondary features of HyS and He I 5876A at high veloc- 
ity (~ 20000 km s"'), produced by complicated and unexplained 
non-LTE effects. Alternatively, using the model atmosphere code 



CMFGE N jHillier & MilleJI 19981) and a ssuming a relevant N en- 
richment. lDessart & Hillien j2005il2006l> reproduce the lines in the 
blue wings of both of Hj8 and He I 5876A detected in the spectra of 
SN 1999em as N II. 

For SN 2005cs a scenario where two line forming regions ex- 
ist for hydrogen and helium is not supported by the line velocities 
measured for the putative HV components. As shown in Fig. |5] 
these line velocities are inconsistent. In particular, the velocity of 
the putative HV \if} component is much larger than that of the HV 
Ho- component. Moreover, the putative HV H/? and HV He I fea- 
tures disappear simultaneously in the spectrum at ~8 days. This 
supports, at least in the case of SN 2005cs, the identification of 
these two features as both due to N II. 

For the same reason, and because of its persistence over 




4500 5000 5500 6000 6500 7000 
rest wavelength 

Figure S. Evolution of the region between 4100 and 7300A. where the min- 
ima of H/?, He I 5876A and Hff are marked, as well as the position of the 
features at ~ 4580A, 5580A and 6300A. 
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Table 3. Model parameters for synthetic spectra shown in Fig. 171 





t= 17'' 


t = 34'' 


log(L/Z.G) 


7.90 


7.88 


Vph [kms"'] 


3710 


2580 


log(Rpi,/Ro) 


3.894 


4.037 


Tph [K] 


7599 


6223 



ious absorption spectral lines in SN 2005cs, as derived from the 
position of their minima. Ha shows the highest velocities, about 
1000 km s"' larger than those of the He I line. The N II and the 
Si II features appear to have smaller velocities than other lines. 
In particular, the velocity of the N II lines decreases from about 
5600 km s"' to ~3200 km s"' between ~3 and 5 days. The veloc- 
ity of the Si II 6355A feature decreases from 4300 km s"' at ~8 
days to about 2100 km s"' one month after the explosion. Over 
the same time interval, the Ho- velocity decreases from 6200 to 
3700 km s"', while the Fe II velocity is ~800 km s"' faster than 
that of the Si II line. In Fig.|51bottom) we compare the Ho- veloc- 
ity ev olution in SN 2005cs with tha t of SN 1987A IPhillips et alJ 
11988), SN 1999em l'P astorellol2o'ol and SN 1999br l Hamuvl200lF 
£astorello et al. 2004J. The Hor velocity curve of SN 2005cs ap- 
pears to be strikingly similar to that of the low- velocity SN 1999br. 
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Figure 6. Top: expansion velocities for Hff, H/3, He I /15876, Na ID, N 
II ^5680, Fe II /i5169. Si II .16355 deduced from the minima of P-Cygni 
profiles in the spectra of SN 2005cs. Bottom: comparison of the Ha velocity 
evolution for SN 2005cs, SN 1987A, SN 1999em and SN 1999br. See the 
text for references. 



a much longer time than the feature at 4580A, we believe that 
the feature near 6300A is Si II 6355A, rather than HV Ho-. The 
presence o f this absorption in the spectra of SN 1999em was 



claimed by'Dessart & Hillieil iim^ . and not explicitly mentioned 
bylSaron et al. ( 200^! 

In Fig.|6jtop) we show the evolution of the velocity of var- 



3.3 Spectral models: the 17 and 34 days spectra 

Some preliminary spectral models have been computed in order to 
provide basic line identification. The code used for the synthetic 
spectra was described i n more deta il by [A bbott & Lucv ( 198fj); 
iMazzah & Lucl <1993h : lLuc\j jl999l) : lMazzalL 1.2000) . The proce- 
dure involves a Monte Carlo (MC) simulation of the line transfer 
based on the Sobolev approximation. The code assumes that all ra- 
diative energy is emitted below a sharp lower boundary. The prop- 
agation of all energy packets is followed through the spherically 
symmetric envelope. Processes of interaction for photons taken into 
account are electron scattering and line transitions. When a photon 
packet is absorbed by a line transition, it is reemitted at a new fre- 
quency corresponding to the branching probabilities for the radia- 
tive decays of the excited level. At the end of the MC calculation a 
formal integral routine derives the emerg ent spectr um based on the 
MC estimate of the source functions ("see lLucvll9 99). 

All m odels were computed assuming solar composition 
iGrevesse & Sauva l 1998). As a starting point for this preliminary 
analysis we used the density structure of the hydrodynam ic ex- 
plosion model adopted to fit SN 1997D fruratto et alll998l) . This 
model was characterised by an outer density profile which can be 
approximated by a power law p oc r"'" at velocities above 3000 
km s"', and contained less than ~ 0.1 Mq of material above this 
velocity. However, in order to reproduce the relatively narrow ab- 
sorptions, in particular the well separated lines of the Ca II IR 
triplet, the density structure was cut above v = 6500 km s"' by 



adopting an even steeper power law (p 



) beyond that ve- 



locity. This modification leads to an insignificant reduction of the 
ejected mass. Tab. |3] gives an overview on the model parameters 
used to derive the synthetic spectra shown in Fig. Input param- 
eters include the total luminosity of the SN and the position of the 
photosphere in velocity space Vph. Together with the epoch t this 
constrains the photospheric radius R and, therefore, the tempera- 
ture at the photosphere via T'^^ 



L/4no-R^ - L/Ancrfv^^^, where 



cr is the Stefan-Boltzmann constant. Tab.|3|gives the final temper- 
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Figure 7. Comparison between synthetic spectra and the observed spectra of SN 2005cs at phase 17 (top) and 34 days (bottom). The parameters used to 
compute the synthetic spectra are reported in Tab.l3l 



ature at the photosphere which is determined iteratively to obtain 
the required emergent luminosity at the outer radius. 

The synthetic spectra in Fig.0are compared with the observed 
spectra at phases ~17 (top) and 34 days (bottom), respectively. The 
bottom panel of Fig.|7|also shows the identifications of prominent 
line features in this spectrum. The synthetic spectra support the 
identification of the feature near 6300 A with the Si II 6347, 6371 A 
lines. In the region between 5000 and 6400 A, together with strong 
Fe II features, prominent lines of Sc II, Ti II, and the Na ID lines 
are identified. 



In the later spectrum we note the increased strength of the 
Ba II lines compared to earlier epochs. Unfortunately, the 5854A 
line is blended with Na ID and the 6497A line with Hor. However, 
the feature at 6I42A is relatively unblended and is detected un- 
ambiguously. The Ba II lines, as well as those of other s-process 
elements (Sc, Sr), are particularly prominent in late-photospheric 
spectra of low-velocity SNe IIP (Pastorello et al. 200 4j) and of SN 
1987 A, while are weaker in other SN e IIP iM azzali et alJll992t 
iMazzali & Chugall995tlutrobin & Ch ueai 200^! 

The proposed density structure does not seem to be able to 



SN2005csinM51 9 




I I I I I I I I I I I I I 

4000 6000 BOOO 10' 

rest wavelength (A) 



Figure 8. Comparison among spectra of SN 2002gd, SN 2005cs, SN 
1999em and SN 1987A shortly after explosion. For references, see the text. 



fit the observations in detail. In particular, it is apparent that the 
proposed density cut, which is required to fit the Ca II IR lines at 
r = 17 d, does not provide a good fit for the epoch / = 34 d. This 
suggests that a steeper density law with less mass at high velocity 
is required to reproduce the low-velocity narrow absorptions seen 
in the observed spectra. 

Concerning the total ejected mass, no conclusion can be drawn 
on the basis of spectral models of early epochs alone. A more de- 
tailed analysis of the density structure would also require informa- 
tion on the later spectral evolution as well as a stud y of the light 
curve, including the duration of the plateau phase (e.g. lTuratto et alJ 
[l998i) . This will be discussed elsewhere. 



3.4 Comparison with other SN IIP 

In Fig.|8|we compare the spectrum of SN 2005cs at ~4 days with 
those of other young SNe IIP: the low -velocity, ^^Ni-poor SN 
2002gd lPastorellQ 2003). SN 1999em iHamuv et al.t2001) and SN 
1987A (Padova-Asiago SN archive). All spectra have blue con- 
tinua and show He I 5876A and the H Balmer lines. Moreover, 
the features we attribute to N II 4623A and 5679A seem common 
to all SN II spectra except SN 1987 A, which has however broader 
lines. 

At ~10 days after explosion (see Fig. |5J, we note a num- 
ber of significant differences between the spectra of SNe IIP. 
While the spectrum of SN 1987A shows a red continuum domi- 
nated by prominent H and Fe II P-Cygni features, the spectra of 
more typical SNe IIP are still relatively blue. Although the spec- 
trum of SN 2005cs is slightly redder than that of SN 1999em 

J HI 

jElmhamdi et a l. 2003), it is bluer than that of the low-velocity, 
^*Ni-poor SN 2003Z I KnoD et al. in oreoaration). The N II features 
visible in Fig.|8|have now completely disappeared, but the promi- 
nent absorption feature at 6300A is now visible in the spectra of 
SNe 2003Z, 2005cs and, as a bump in the Ha abso rption profile, 
in 1999em (see discussion in lElmhamdi etai]|2003^ . The identifi- 
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Figure 9. Comparison among spectra of SN 2003Z, SN 2005cs, SN 1999em 
and SN 1 987 A at a phase of ~ 1 days. A blow-up of the Ho' region is shown 
in the top-right corner, showing the different strengths of the Si II feature 
in SN IIP spectra. 
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Figure 10. Spectra of the low-velocity SN IIP 1999br, 2003Z and 2005cs 
at a phase of about two weeks. Spectra of SN 1999em at a phase of 17 and 
21 days are included for comparison. Again, a blow-up of the Ho- region is 
also shown in the top-right comer. 

cation of this feature as Si II 6355A is supported by the spectral 
modelling presented in Sect. l3.3l 

In Fig. ^1 we compare the spectrum of SN 20 05cs at 
phase ~ 2 weeks to those of the ^^Ni-po or SNe 199 9br (Hamuvj 
2001 ; Pastorello et al. 2004) and 2003Z ( Kn op et aljin oreDaratio J 
at a similar epoch. We also include 2 spectra of SN 1999em 
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Figure 11. Evolution of the continuum temperature of SNe 2005cs, 1987A, 
1999em, 1999br. For references, see text. 

jLeonard et allbOOl iHamuv et al]l200ll) . taken about 17 and 21 
days after the explosion. The Si II 6355A feature is prominent in 
all spectra, included that of SN 1999em at phase ~I7 days, but it is 
not longer visible 4 days later. 



4 DISCUSSION 

The observationa l analogies between SN 2005cs and objects simi- 
lar to SN 1997D <Turattoe tal.°1998') are remarkable. In Tab.|3|we 
report some significant information available in the literature for a 
number of objects belonging to this group. Only SNe with very low 
ejecta velocities and/or small ejected ^*Ni mass 10"^ Mq) are in- 
cluded in Tab. |4| All these objects belon g to the faint tail of the 
luminosity distribution of SNe IIP ( Ham uvl200 ll: |Pastorellol2003t 
lHamuvl2003tlPastorello et ali2005i : iZampieril2005h and have sim- 
ilarly low expansion velocities. In particular, the good match be- 
tween the Hff velocity curves of SNe 2005cs and 1999br (Fig.|6| 
bottom), provides further support to our idea that SN 2005cs can 
be regarded as another S N 1997D-like event. In an alogy to other 
low-luminosity SNe IIP jPastorello et aljlin preoJ) . we therefore 
expect that SN 2005cs will also evolve through a long plateau (3-4 
months) and reach a faint late-time luminosity. This would be in- 
dicative that a very small mass of '^Ni (of the order of ~ IO^^Mq or 
less) was ejected in the explosion, as the other similar SNe of Tab. 

H 

The earliest spectra of SN 2005cs are very blue, suggesting a 
very high continuum temperature (2-3x10'* K). In Fig. ll li the evo- 
lution of the temperature of SN 2005cs derived from a blackbody 
fit to the spectral co ntinuum is compared with th ose of SN 1987 A 
f Phillips etal.'l988"), SN 1999em (as measured in'Pastorello'2003') 
and SN 1999br (Pastorello et al. 2004). Again, the continuum tem- 
perature evolution of SN 2005cs is not different from that of other 
SNe IIP. In particular, SN 2005cs has higher continuum tempera- 
tures than SN 1999br, but slightly lower than SN 1999em. Inde- 
pendent on the phase, all these long-plateau SNe have significantly 



higher continuum temperatures than SN 1987 A. The temperature 
of SN 1987 A becomes stationary about 3 weeks after the explo- 
sion, while for other SNe IIP the same happens 30-40 days past 
explosion. This indicates that in SN 2005cs (and in the other SNe 
IIP) the H envelope starts to recombine later than in SN 1987 A. 

If the faint absolute luminosity and the very well developed 
plateau should be confirmed by the observational campaign still 
in progress, they would fit with some difficulty the relatively small 
mass QAzAMS ~ 9Mh) progenitor scenario derive d from the progen- 
itor detection iMaund et all2005l : lLi et all200'3) . As mentioned in 
Sect.0 both groups argued that the low luminosity of the progeni- 
tor indicates that it had a small mass. Maund et al. ( 2005) obtained 
a likely range of bolometric luminosity 4 < log(L/LQ) < 4.4 (shaded 
region of their Fig. 3) , and even for conservative errors, log(L/Lo) 
< 4.6. iLi et alJ i2006h obtained an even lower value. To compare 
with the observed luminosity, they used the Geneva evolutionary 
models. The highest luminosities attained in the Geneva models 
are log(L/Lo) = 4.2, 4.5, and 4.8 for M = 7Mo, 9Mo, and 12Mq, 
respectively. To be consistent with the observed luminosity for the 
conservative limit, they suggested M = 7-12Mo for the progenitor. 

However, the Geneva evolutionary models do not reach the 
pre-supernova stage. They cover the evolution up to the formation 
of 0+Ne+Mg core for M > 9Mo and only up to the formation of 
the C+O core for M = 7Mq. Stars with M < SMq form a degenerate 
C-l-O core whose mass yicore increases t owards the Chan drasekhar 
mass, Mrfi, as the star climbs the AGB ('PaczvnskT 1970). Eventu- 
ally the stars will either lose their H-rich envelope to form C-l-O 
white dwarfs or undergo a thermonuclear explosion (the so-called 
Type 1-1/2 supernovae) when M„;e gets close to Mrf, (for reviews, 
see Sugimoto & Nomoto 1980; Nomoto & Hashimoto 1988). The 
8-lOMo stars form a degenerate 0-l-Ne+Mg core whose M.core also 
increases toward Mrf, on the AGB liNomoto«1984i') . Their final fate 
is also either an O-l-Ne-l-Mg white dwarf or a core-colla pse super- 
nova when Mcorf gets close to M^fc jNomotJl98ill987h . 

The luminosity L of AGB stars withM < lOM© obeys Paczyn- 
ski's (1970) Mco„-L relation. If a star with M < IOMq reaches 
M-core = 1 4Mq, the pre-supernova luminosities should be at least 
log(L/Lo) = 4.8 for X(H) = 0.7, and e ven higher if the He ab un- 
dance is higher (e.g. because of mixing. iHashimoto et al Il99l). If 
the progenitor was as massive as 12Mq, the pre-supernova lumi- 
nosities would be log(L/Lo) > 4.8. Systematic studies of the pro- 
genitor luminosity, including the effect of rotation, would be desir- 
able to obtain a more reliable detailed comparison. 

Therefore, the observed luminosity of the putative progenitor 
of SN 2005cs is inconsistent with the pre-supernova luminosity of 
7-12Mq stars. If large extinction causes the observed luminosity of 
the progenitor to be very small, the observed luminosity cannot be 
used to constrain the progenitor's mass. One possibility is that the 
absolute magnitude (and hence the mass) of the progenitor star is 
underestimated because of dust enshrouding the supergiant progen- 
itor a nd later swept away by the SN explosion fcl raham & MeikI3 
Il98d) . This scenario was ruled out bv lMaund et alJt200 5) because 
no K band excess in the magnitude of the progenitor was detected 
in archival images. However, we cannot exclude the possibility that 
a particular dust composition and an extremely low temperature 
cause the light to be absorbed at optical waveleng ths and re-emitte d 
mostly in the mid- and far-IR bands (see e.g. Pozzo et al.'"2004|). 
Dust was ob served in the la te evolution of another well studied SN 
up 2003gd faendrv et al] |2004V Despite having a plateau lumi- 
nosity and an expansion velocity evolution similar to the "normal" 
SN 1999em, SN 2003gd showed low late time luminosity, indicat- 
ing the ejection of a small mass of ^''Ni (0.015 Mo), only a factor 
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Table 4. Main data about low-luminosity, Ni-poor SNe IIP available in literature. In column 5 we reoprt the expansion velocity of the ejecta at phase ~35 
days, obtained from the position of the Ha P-Cygni absorption. 



SN name 


Host Galaxy 


Galaxy Type® 


Mv,,,* 


vg' (kms-') 


^i^Ni Mass (Mo)* 


Source 


1994N 


UGC 5695 


Sab 


-14.8 


4250 


5x10-3 


♦ 


1997D 


NGC 1536 


SBc 


-14.3"^ 




7x10-3 


♦ 


1999br 


NGC 4900 


SBc 


-13.5 


3700 


2x10-3 


♦ 


1999eu 


NGC 1097 


sBb 


-13.7"^ 




<3xlO-3 


♦ 


2000em 


LEDA 143614 


7 


-16.3* 




unknown 


★ 


200 Idc 


NGC 5777 


Sb 


-14.1 




5x10-3 


♦ 


2002gd 


NGC 7537 


Sbc 


-15.6 


4100 


<3xlO-3 


■ 


2003Z 


NGC 2742 


Sc 


-14.6 


4000 


5x10-3 


■ 


2004cm 


NGC 5486 


Sm 


-13.9* 




unknown 


* 


2004eg 


UGC 3053 


Sc 


-14.7:* 




unknown 


▲ 


2005cs 


M51 


Sbc 


-15.2 


3700 


unknown 


T 



♦iPastorello et al. 1 2004); ★ Strolaer et al. 1 2000); ■ Pastorello 1 2003); A Youns et al. 1 2004); FiliriDenko et al. 1 2004); T this paper. 
* http:iicheopslMchicago.edu/pub/mefchart/meCand-53088-1324-491-run003712-20-5-0189-00031.html see also .ConnoUv et alj|2004l) . 
® LEDA; ' Absolute magnitude measured at the end of plateau; * Unknown phase (but during plateau); R band absolute magnitudes. 
* fi computed assuming Ho = 72 km s-' Mpc-' ; when possible. My,,/ was computed about 1 month after the explosion. 



two m ore than SN 1997D. IVan Dvk et alj iiooal) and lSmartt etall 
<2004l) found a moderate-mass progenitor for SN 2003gd (8 Mq). 
However, an evident light echo was recently detected in HST im- 
ages of this object ( Van Dvk et al. 2006; Suaerman 2005). This was 
due to SN light scattered by large-grain dust, located 110-180 pc 
in front of the SN, that survived the initial UV-flash. If this mate- 
rial was present before the SN explosion, it could be responsible for 
significant extinction of the star light, leading to an underestimate 
of the luminosity (and hence of the mass) of the progenitor. 

As further support to this discussion, most SNe reported in 
Tab. 13 exploded in late Type galaxies (mainly of morphological 
Type Sb-Sc, by LEDA'), which would be consistent with massive 
progenitors. 

Finally, there is some ambiguity as regards the location of the 
putative progenitor of SN 2005cs. Therefore, even if the works of 
[Maund et al. ( 2005) and Li et al. ( 2006) seem to support a moderate 
mass scenario for the precursor star of SN 2005cs, late time obser- 
vations of the explosion site will probably be required in order to 
remove the residual uncertainty in the correct identification of the 
progenitor. 
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